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Abstract-Ill  is  paper  presents  two  methods  to  estimate  the 
two  dimensional  (2-D)  direction  of  arrival  (DOA)  for 
coherent  and  non-coherent  sources.  The  proposed  methods 
have  many  advantages  over  existing  schemes.  First,  they 
construct  the  data  from  a  single  snapshot  in  a  Toeplitz  form, 
whose  rank  is  directly  related  to  the  DOA  of  signals, 
whether  the  signals  are  coherent  or  not;  hence,  the 
algorithm  does  not  require  any  forward/backward  spatial 
smoothing.  Second,  the  two  proposed  methods  can  rapidly 
estimate  the  2-D  DO  As  of  incident  signals  without  requiring 
singular  value  decomposition  (SVD)  or  eigenvalue 
decomposition  (EVD),  even  in  the  case  of  coherent  signals 
and  a  single  snapshot.  The  two  methods  are:  (1)  orthogonal 
projection  real-valued  rank  revealing  QR  factorization  (OP- 
RRRQR),  and  (2)  orthogonal  projection  real-valued  rank 
revealing  LU  factorization  (OP-RRRLU).  The  proposed 
methods  reduce  computational  complexity  and  the  cost  at 
least  by  a  factor  of  four  by  applying  a  unitary 
transformation,  to  the  complex  Toeplitz  form  to  real  data 
without  forming  the  covariance  matrix.  The  proposed 
algorithms  employ  the  unitary  root  MUSIC  and  unitary 
MUSIC  using  cross  array  configuration  to  estimate  the  2-D 
DOA  azimuth  and  elevation  angles  without  using  the 
extensive  2-D  MUSIC  search.  Hence,  they  can  reduce  the 
computational  load  and  cost  significantly  and  can  be  applied 
in  real-time  radar/sonar  and  commercial  wireless  systems. 
The  simulation  results  show  that  the  proposed  algorithms 
can  efficiently  estimate  the  2-D  DOAs  from  different  sources. 

i  Tli is  work  was  supported  by  a  research  contract  from  Air  Force  Research 
Laboratory/Clarkson  Aerospace  under  Prime  Contract  Number  FA865 0-05 -D- 
1912. 

I.  INTRODUCTION 

The  problem  of  estimating  the  two-dimensional  (2-D) 
direction  of  arrival  (DOA),  azimuth  and  elevation  angles, 
for  incident  signals  on  an  antenna  array  has  received 
considerable  attention  in  the  field  of  array  signal 
processing  [1-13].  This  problem  has  applications  in  many 
fields  including  radar,  sonar,  radio  astronomy,  seismic  data 
processing,  and  mobile  communication  systems.  Existing 
algorithms  employ  either  eigenvalue  decomposition  (EVD) 
of  the  sample  covariance  matrix  or  singular  value 
decomposition  (SVD)  of  the  data  matrix  to  estimate  the 
signal  and  noise  space,  and  are  computationally  extensive 
and  time  consuming,  especially  when  the  number  of 
antenna  array  elements  N  is  larger  than  the  number  of 


incident  signals.  Furthermore,  many  of  these  methods 
require  a  two  dimensional  search,  complex  pair-matching 
of  the  azimuth  and  elevation  angle,  have  an  estimation 
failure  problem  when  the  elevation  angles  are  between  70° 
and  90°,  need  a  large  number  of  samples  to  estimate  the 
covariance  matrix,  and  fail  when  the  sources  are  highly 
correlated  or  completely  coherent. 

Die  2-D  DOA  estimation  problem  in  [1]  requires  an 
exhaustive  2-D  peak  search  through  all  possible  steering 
vectors.  In  addition,  the  sources  must  be  non-coherent  and 
a  large  number  of  snapshots  are  required  to  estimate  the 
covariance  matrix.  In  [2-5]  the  authors  proposed  a  2-D 
DOA  estimation  method  using  a  different  array 
configuration.  While  this  method  does  not  require  a  2-D 
peak  search,  it  has  other  drawbacks:  (1)  it  requires  a  pair- 
matching  between  the  2-D  azimuth  and  elevation  angle 
estimations,  and  (2)  estimation  fails  when  the  elevation 
angle  is  between  70°  and  90°  or  when  the  signal-to-noise 
ratio  (SNR)  is  Low. 

Die  methods  in  [1-4,  6-13]  require  a  large  number  of 
snapshots  to  estimate  the  covariance  matrix.  Diis 
covariance  matrix  will  have  a  Toeplitz  structure  if  the 
incident  sources  are  uncorrelated  and  statistically 
stationary.  However,  if  the  incident  sources  are  coherent 
(e.g.,  in  a  multipath  environment),  the  covariance  matrix  is 
not  Toeplitz.  In  order  to  obtain  the  Toqilitz  structure,  the 
preprocessing  spatial  smoothing  technique  [14-16]  has 
been  introduced  in  order  to  destroy  the  cross  correlation 
between  directional  components;  however,  smoothing 
technique  requiresa  large  number  of  snapshots  and 
averaging  the  covariance  matrices  which  increase  the 
computational  load  and  the  complexity. 

In  this  paper,  we  present  a  new  approach  to  solve  the 
2-D  DOA  estimation  problem.  There  are  five  main 
underlying  ideas  behind  this  approach.  First,  the  proposed 
algorithm  preserves  the  Toeplitz  structure  by  mapping  die 
(2V+l)xl  data  vector  from  a  single  snapshot  into  an 
(N+l)  x  (N+l)  Toeplitz  matrix  whose  rank  is  related  to  the 
DOA  of  the  incoming  signals  independent  of  whether  the 
sources  are  coherent  or  not.  Hence,  die  proposed  method 
does  not  use  forward/backward  spatial  smoodiing  resulting 
in  a  reduction  in  the  computational  complexity  and  cost. 
Second,  the  proposed  method  uses  only  a  single  snapshot 
of  the  received  signals  to  estimate  the  DOA  of  the  incident 
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sources.  This  reduces  the  computational  load  drastically 
and  makes  the  proposed  method  a  good  candidate  for  real¬ 
time  implementation.  Third,  to  reduce  the  computational 
complexity  and  cost  further,  we  employ  the  unitary 
transformation  in  [17]  to  convert  the  complex-valued 
Toeplitz  data  to  real-valued  data.  This  reduces  the 
processing  time  by  almost  a  factor  of  four,  since  the  cost  of 
complex  multiplication  is  four  times  that  of  real 
multiplication.  Fourth,  the  two  proposed  methods  can 
rapidly  estimate  the  2-D  DO  As  of  incident  signals  without 
requiring  SVD  or  EVD,  even  in  the  case  of  coherent 
signals  and  a  single  snapshot.  The  two  methods  are:  (1) 
orthogonal  projection  real-valued  rank  revealing  QR 
factorization  (OP-RRRQR)  and  (2)  orthogonal  projection 
real-valued  rank  revealing  LU  factorization  (OP-RRRLU). 
The  rank  revealing  QR  factorization  [18-20]  is  precisely 
used  to  compute  the  subspace  information  and  effectively 
update  the  signal  information  that  can  be  used  to  track  a 
moving  source.  The  QR  factorization  is  also  applied  to 
estimate  the  DO  As  in  the  non-stationary  environment  of 
tracking  moving  sources.  Fifth,  we  propose  a  cross  array 
configuration  that  consists  of  two  centro-symmetric 
uniform  linear  arrays  in  the  z-x  plane;  then,  the  proposed 
unitary  root  MUSIC  is  applied  to  the  array  in  the  z  axis  to 
estimate  the  elevation  angle.  Subsequently,  the  unitary 
MUSIC  method  is  used  to  estimate  the  azimuth  by  the  use 
of  the  estimated  elevation  angle  for  each  source.  Finally, 
the  proposed  methods  do  not  require  any  pair-matching 
between  azimuth  and  elevation,  they  do  not  suffer  from 
estimation  failure,  and  avoid  the  2-D  search  peak. 

The  rest  of  this  paper  is  organized  as  follows.  Section 
II  presents  the  proposed  2-D  DOA  methods.  Section  III 
shows  simulation  results,  and  in  Section  IV  we  draw  our 
conclusions. 


II.  SYSTEM  MODEL 

A-Orthogonal projection  with  reed-valued  rank-revealing 
QR  factorization 

Figure  1  shows  the  proposed  array  configuration 
consisting  of  two  symmetric  uniform  linear  arrays  (ULA) 
with  interspacing  d  equal  to  a  half  wavelength  of  incident 
signals.  We  assume  that  all  the  incident  sources  have  the 
same  carrier  frequency.  Each  uniform  linear  array  in  Figure 
1  consists  of  2N  elements,  and  the  element  placed  at  the 
origin  is  numbered  0  for  reference  purposes.  One  array  is 
placed  on  the  z  axis  and  the  other  on  the  x  axis.  Suppose 
that  there  are  K  narrow  band  sources  where  the  A'-th  source 
has  an  elevation  angle  0k  and  an  azimuth  angles ,  k=l, ... , 
K. 

Step  1.  Estimation  of  Elevation  Angle  0k 

For  a  given  snapshot  t ,  the  output  signal  from  the  k- th 

element  on  the  z  axis  is  given  by 


mo- +-.»(<)  id 

i=0 

where  si  (t )  is  the  signal  from  the  /-th  incident  source,  and 
rifr  (?)  is  the  noise  at  the  A-th  element. 

If  we  use  the  element  at  the  center  of  the  array  as  a 
reference  point,  then  the  {2N  + 1)  x  1  output  vector  from 
the  2N  + 1  antenna  elements  placed  on  the  z  axis  can  be 
written  as 


z(,)= 


=  A(0)  S(t)+N(t) 


_ZA() 

where 


A(^)  =  [a(^)  a(02)  ...  a<0r)] 
is  the  i  2N  +  l)x  K  array  response  matrix, 
where 


(2) 


(3) 


a(ek)= 

kT  -  ■  -  < 

T 

(4) 

is  die  corresponding  (2V  +  l)x  1  array  response  vector, 

with 

uk  =  exp  {-jin  cl  cos  ( 0k )/  X) , 

(5) 

where  S(t)  is  the  vector  of  received  signals 

S«)  =  [i 

htt)  si(0  MOf 

(6) 

and 

tf(0  =  [ 

(«_Ar(0  «o(0  —  - 

(7) 

where  N(t)  is  i 

the  (2N  + 1)  x  1  noise  vector.  Herein, 

die 

superscripts  T  and  *  denote  the  transpose  and  conjugate 
operations,  respectively. 

In  the  proposed  method,  we  map  the  output  data  vector 
Z(t)  with  dimension  (2Ar  +  l)xl  to  a  Toeplitz 
Hennition  data  matrix  with  dimension  (N +l)x  (N+ 1)  . 
The  advantage  of  introducing  the  Toeplitz  Hennition  data 
matrix  is  that  it  has  a  rank  that  is  related  to  the  DOA  of  the 
sources  independent  of  whether  the  sources  are  coherent  or 
not.  This  matrix,  denoted  Y,  is  given  by 


% 

2-1 

2-2 

Z-N 

zi 

20 

2-1 

Z2 

Z1 

20 

Z-(W-2) 

Jn 

ZN- 1 

ZN-l 

20 
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where  we  have  dropped  the  time  index  t 

Assuming  a  noise  free  case  the  Hermitian  Toplitz  data 

matrix  in  (8)  can  be  rewritten  in  terms  of  7?(6>),  ,  and 

S  as  follows 

y=|b(^s  b(<9)<d;(0)s 

where 

B(«9)=[bC6o  h&2)  ...  b(eK)] 

is  the  ((AT  + 1 )  x  K)  array  response  matrix  where 


Y  =  B{0)DB{0)H  (16) 

Unitary  transformation  Method  for  Toeplitz  data 
Let  J  represent  the  exchange  matrix  (i.e.,  l’s  on  the 
antidiagonal  and  0’s  elsewhere)  as 


B( 


to 

i _ i 

'vo 

J  = 

r0 

vl 

(10) 


"WH1 

and 


0 

U2 

0 

0 


-  diag(i fj* 


(11) 


uK)  (12) 


The  square  Hermitian  Toeplitz  in  (10)  can  be  decomposed 
as  follows 


Y  =  WDW 


H 


where 


W  = 


and 


D  = 


(13) 


'  1 

1 

1 

1  " 

U3 

■  UK 

“ l 

i4 

(14) 

U2 

U3 

-  < 

0 


0 

0 

0 

sz 

rH 


(15) 


Since  7  -WDWn  and  the  diagonal  signal  matrix  D 
is  of  full  rank,  the  rank  of  7  is  the  same  as  that  of  W .  The 
matrix  W  in  (14)  has  the  structure  of  an  (N  + 1)  x  K 
Vandermonde  matrix.  The  rank  of  W  is  equal  to  the 
Minimum  of  (N  + 1)  and  K ;  hence,  rank  of  W  is  equal 


(17) 


Note  that  J  =  JH  =J~l 

Toeplitz  matrices  are  persymmetric ,  meaningthat  they  are 
symmetric  about  their  southwest-northeast  diagonal.  For 
such  a  matrix  P 

jptj  =  p 

(18) 

Let  U  be  an  MXM  matrix  defined  as 


G  = 


1 


i  j 
jJ  -fi 


(19) 


G  =  ~T 
V2 


where  J  and  /  are  exchange  and  identity  matrices, 
respectively,  with  dimensions  (M/2)x  (M/2),  and  G  is  a 
unitary  matrix  which  satisfies 

G* J  =  G .  (20) 

Equation  (19)  holds  for  an  M  even.  In  the  case  where  M  is 
odd,  the  unitary  matrix  G  can  be  written  as 
"7  0  J 

0  Yl  0  (21) 

JJ  0  ~JI. 

where  /  and  J  have  dimensions  of 

((M-l)/2)x((M-l)/2)  and  0  =  [0,0,...,0]  . 

Construct  the  output  data  of  the  proposed  method  as 
follows 

+  (22) 

If  we  pre-multiply  the  construct  data  in  (22)  by  G ,  post- 
multiply  by  Gl!  and  take  die  conjugate,  we  get 

(G'PGH)*  =  [G(7  +  J7V)GHJ 

=  (G(7+7ff)Gff)*  =G*7*Gr  +  GTrG7'  (23) 


to  K .  This  means  that  die  rank  of  the  Hermition  Toeplitz 
data  matrix  is  equal  to  the  number  of  DO  As  of  the  sources 
whedier  the  sources  are  coherent  or  not.  Therefore,  all  the 
incident  sources  can  be  detected  even  if  the  sources  are  not 
coherent  because  die  Toeplitz  data  matrix  structure  will  be 
preserved  in  bodi  scenarios  (i.e.,  coherent  and  non-coherent 
sources). 

If  W  matrix  in  (13)  is  the  same  as  die  collection  of 
array  response  vectors  from  different  directions  in  B{o)  in 
(11),  then  equation  (13)  can  be  rewritten  as 


*  j-j 

Note  that  because  7  is  a  Toeplitz  matrix  dien  JY  J  =7  . 

Using  JJ=I ,  we  can  rewrite  (23)  as  follows 

(GxVIJh)  =  G'JJY'JJG7  +  G'JJYtJJGt 

=  G*J(JY  *J  )JGt  +  G'J  [JY  TJ  )JGt 

=  G’J(Yh  +  Y)JGt  (24) 

If  we  substitute  (20)  for  (24),  we  get 

(G'FG")*  =g(7  +  7")g* 

=Gx¥Gh  (25) 
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Therefore,  Gy¥GH  is  a  real-valued  matrix;  the 
decomposition  of  GT'G^  requires  only  a  real  computation 
which  means  the  computational  load  and  cost  will  reduce 
significantly  without  effecting  the  accuracy  of  the  DO  As. 

The  proposed  method  uses  the  orthogonal  projection 
real-valued  rank  revealing  QR  factorization  to  estimate  the 
2-D  DO  A  elevation  and  azimuth  angle  from  coherent  /non¬ 
coherent  sources  by  using  the  real  data  matrix  in  (25).  The 
necessary  information  about  the  noise  subspace  or  the 
signal  subspace  can  be  extracted  using  the  rank  revealing 
QR  factorization  [18-20].  One  of  the  reasons  that  QR 
factorization  is  widely  used  in  adaptive  applications  is  that 
in  RRQR  the  signal  information  can  be  effectively  updated, 
making  the  algorithm  suitable  to  track  moving  sources. 

The  real  data  matrix  G'f'G"  in  (25)  with  dimensions 
(N  +  l)x  (iV  +  l)  using  the  QR  factorization  can  be 
expressed  as  the  product  of  a  real  and  orthogonal 
(N  + 1)  x  (N  + 1)  matrix  and  an  (N  + 1)  x  (N  + 1)  rank- 
revealing  upper  triangular  matrix  with  real  entries.  Then 


the  factorization  of  G'i'G^  can  be  written  as 


Gx¥Gh  =  QR  =  Q 


R, 


R 

R 


12 

22 


where  Rn  is  an  ( KxK )  matrix, 
an  K  x  (N — K  + 1)  matrix  and  R12 


(26) 

Rn  is 
is  an 


(N-K  +l)x(N-K  + 1) matrix,  and  all  real  matrices.  The 


QR  factorization  in  (26)  is  called  rank-revealing  QR 
factorization  if  R22  has  a  small  norm.  Since  R22  is  very 
small,  the  basis  of  the  noise  space  can  be  obtained  from  the 
above  R  factor.  Let  V  represent  any  vector  in  the  null 
space  of  R ,  i.e.,  RV  =  0.  To  find  the  structure  of  V  9  we 
partition  V  into  v1  with  K  components  and  v2  with  (N- 
K+l)  components.  Then  RV  -  0  implies  that 


[*n  *n]  1 


=  0 


(27) 


so  that  Rnv,  +  R12v2  =  0  .  Since  R^  is  a  non-singular 


matrix,  can  be  written  in  terms  of  v2  as  follows 


vi  =  ~RnRuv2 

Thus  V  can  be  written  as 


(28) 


v  = 

V 

3. 

_  I  N-K+l 

(29) 


To  find  the  basis  of  the  null  space  of  R,  we  choose  any 
set  of  N  -K  + 1  linearly  independent  vectors;  for  example, 
the  columns  of  the  N  -  K  + 1  dimensional  identity  matrix. 
The  basis  for  the  null  space  of  the  upper  triangular  matrix 


R,  which  is  also  the  null  space  ofG'f'G^ ,  is  therefore  given 
by 


E  = 


~^Rn 


(30) 


It  is  important  to  observe  that  the  columns  of  Ez  are 
not  orthonormal  in  contrast  to  the  null  space  which  can  be 
derived  from  the  SVD  or  EVD  techniques.  We  now  note 
that  the  subspace  spanned  by  the  columns  of  Ez  is 
orthogonal  to  the  subspace  spanned  by  the  columns 
of  B(O)  ,  where  the  column  of  B{0)  contains  the 


information  about  the  AO  As  of  incident  sources.  This  is 
similar  to  the  well-known  MUSIC  algorithm  in  which  the 
eigenvector  of  the  noise  subspace  is  orthogonal  to  the 
steering  vector  of  the  signal  subspace.  To  find  the  AO  As, 


we  search  the  minimum  peaks  of 


ErB{6) 


.  Since  the 


basis  of  Ez  is  not  orthononnal  we  use  the  orthogonal 
projection  onto  this  subspace  to  improve  the  performance 
by  making  the  basis  of  the  null  space  of  Ez  orthonormal. 

Since  we  consider  die  uniform  linear  array,  die 
proposed  method  employs  the  root  MLTSIC  [21]  to 
estimate  die  AO  As  for  the  incident  signals.  The  proposed 
unitary  root  MUSIC  converts  the  power  spectrum  of  the 
MUSIC  algorithm  into  a  polynomial  whose  roots  contain 
information  about  the  elevation  angle  0  as 

F~l  ( 0)  =  Bh  ( 0)  GhEzEzhGB  (0)  (31) 

Since  the  RRQR  factorization  is  for  die  real- valued 
matrix  instead  of  the  complex-valued  matrix  in 
conventional  Root  MUSIC,  then  the  complexity  of  the 
proposed  mediod  Unitary  Root  MUSIC  in  (31)  is  about 
four  times  lower  than  the  conventional  Root  MUSIC. 

To  make  the  basis  of  Ez  orthononnal  we  use  the 
orthogonal  projection  onto  this  sub  space,  which  is  given 
by 

W0=Ez(e?EzYe»  (32) 

This  implies  that  [14] 

w0B{e)=  0  (33) 

Now  equation  (31),  using  orthogonal  projection,  becomes 
F-;  =  B{6)HGHWQGB(e)  (34) 

The  roots  of  the  polynomial  in  (34)  can  be  used  to 
estimate  the  elevation  angle  0 ^  of  the  incident  signals. 

Step  2.  Estimation  of  Azimuth  Angle  (j)k  with  Estimate  0k 

The  estimate  0 ^  obtained  in  Step  1  will  be  used  for 
estimation  of  die  azimuth  DOA  (j)k 

The  proposed  orthogonal  projection  real-valued  rank 
revealing  QR  factorization  will  employ  the  signal  vector  X 
received  at  the  ULA  in  the  x  axis  direction  and  the  MUSIC 
algorithm.  By  doing  diis  we  can  avoid  the  failure  problem 


4 


Authorized  licensed  use  limited  to:  AFRL  Technical  Library.  Downloaded  on  August  18,2010  at  15:03:47  UTC  from  IEEE  Xplore.  Restrictions  apply. 


in  the  joint  azimuth  and  elevation  angle  estimation  and 
pail -matching  problem. 

The  ((2N+l)x  1)  signal  vector  received  from  the 
symmetric  uniform  linear  array  which  is  a  function  of 

{pi  <j>i )  on  x-axis  is  given  by  ...  Xq,---,xnY  .  In  the 


proposed  method  we  map  the  ( (2N+  l)xl)  vector  to  the 
(N  + 1)  x  (N  4- 1)  Toeplitz  data  matrix  as  follows 


(35) 


*0 

x_x 

x_2 

X-N 

* 

x0 

X-! 

X-(N- 1) 

x  = 

*2 

Xi 

xo 

X-(N- 2)  Z 

XN- 1  • 

XN- 2 

xo 

The 

data 

matrix 

in 

(38) 

can  be  r 

O *9S9  and  C(09<f)  as 


c{{q,<p)s  c{e,cp)^2{e,<p)s  ■■■  c(6,<p)(&2{&,<pjfMs 


+Qx 


C(0,<f>)  =  {c(8h<fa),  c{Q2,<f>2),  •  c(%,^ ))  (36) 

(6k,<pk)=(l,  (ek),  ...,  (, ekf )  (37) 

.  2nd  sin  0k  cos  <pk 


%=exP  ~J 


<D-  = 


X 

e[  0  •••  0 

0  e\  ■■■  0 

0  0 
0  0  -  * •  et 


=  diag(e1m  e*  eK) 


(38) 


where  Qx  is  a  complex  noise  variable.  Using  the  unitary 
transformation  we  convert  the  complex  data  matrix  in  (36) 
to  real  data  as  follows 

T  =  X  +  JXmJ 

GYGh  =  g(x  +  JX’j)gh  (39) 

By  applying  the  same  rank  revealing  QR  factorization 
procedure  in  (26)-(34)  on  GYGff,  the  null  space  Z^can  be 

found  in  the  same  way.  The  orthogonal  projection  onto  this 
subspace  is  given  by 

Q0=Ex[ExExYe?  (4°) 

The  azimuth  angle  estimation  (j>  can  be  obtained  using 


the  estimation  0  which  is  already  found  in  Step  1.  By 
employing  the  proposed  Unitaiy  MUSIC  search  peak  (j>^ 


can  be  found  from  the  maximum  peaks  of  the  following 
power  spectrum  as 


P 


_ 1 _ 

c(ek,<p)H  gh  EX  Gc(ek,<p) 


(41) 


for  source  k :9  k=l,  K  wdiere  Ex  is  similar  to  Ez  in  (30). 
Note  that  the  2-D  DOA  dimensional  search  reduces 
approximately  to  a  1-D  search  since  the  number  of  sources 
is  very  small.  In  addition,  die  complexity  of  the  proposed 
Unitaiy  MUSIC  method  in  (41)  is  about  four  times  lower 
than  die  conventional  MUSIC  [24]. 


B-  Two-Dimensional  DOA  Estimation  with  orthogonal 
projection  and  real  rank  revealing  LU  factorization 

In  this  section  we  present  a  metiiod  to  obtain  the  2-D 
DOA  estimation  azimudi  and  elevation  angles  from 
coherent  non -coherent  sources  using  real  rank  revealing 
LU  factorization.  This  method  is  referred  to  as  orthogonal 
projection  real  rank  revealing  LU  factorization  (OP- 
RRRLU).  Rank  revealing  LU  factorization  [23]  can  reduce 
the  complexity  over  that  of  (RRQR)  by  a  factor  of  two. 

The  RRRLU  factorization  of  GT'G*  can  be  written  as 


U(YY)luUH  =LU  =  L 

where  Un  is  an 
an  Kx{N-K  + 1) 


"^11 

i 

CN 

b  ! 

(42) 

0 

U22_ 

(K  x 

K) 

mahix, 

Ui2  is 

matrix 

and  U2  2 

is  an 

(N -K +  l)x(N -K +1)  matrix,  and  wdth  real  entries. 
Since  Un  is  a  non-singular  matrix  widi  full  rank  equal  to  K 
and  U2 2  has  a  small  norm,  the  factorization  in  (42)  is 
referred  to  as  real  rank  revealing  LU  factorization.  The 
basis  of  the  noise  space  can  be  obtained  from  U  as  follows: 
the  upper  triangular  matrix  U  represents  die  null  space  of  Y. 
Since  U22  has  a  very  small  noun,  the  upper  triangular 
matrix  U  can  be  written  as 

U  =  [Un  Un]vT  (43) 

where  W  is  a  permutation  matrix  that  represents  die  row 
and  column  interchanges.  Let  q  represent  any  vector  in  the 
null  space  of  U,  i.e.? 

Uq=  0  (44) 


To  find  the  structure  of  q,  we  partition  q  into  qx  with  K 


components  and  q2  with  (N-K+l)  components.  Equation 
(43)  can  be  written  as 


Un 


=  0  =>  Unq1+Unq2  =0  (45) 


From  (45)  and  the  fact  that  11^  is  full  rank  with  rank  K, 
qx  can  be  written  as 


ft  =  ~uiiun<i2 

Using  Equations  (47-48)  vector  q  can  be  written  as 


q  = 


XX 


<h 


(46) 

(47) 


To  find  the  basis  for  the  null  space  of  U9  we  choose  any  set 
of  N-K+l  linearly  independent  vectors.  The  most  obvious 
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choice  is  the  set  of  linearly  independent  columns  of  the 
identity  matrix.  A  basis  for  the  null  space  of  the  upper 
triangular  matrix  U9  also  the  null  space  of  7,  is  therefore 
given  by 


Q  = 


-UnUu 


(48) 


Note  that  columns  of  Q  are  not  orthonormal  as  the  ones 
provided  by  SVD  or  EVD.  We  observe  that  the  null  spaces 
of  (OP-RRRQR)  in  (30)  and  the  one  in  (48)  for  the 
proposed  OP-RRRUL  have  the  same  form  and  dimensions. 
We  use  the  orthogonal  projection  onto  this  subspace  to 
improve  the  performance  by  making  the  basis  of  the  null 
space  of  V  orthonormal.  The  columns  of  the  null  space  of 
V  become  orthonmal  as  follows 


Q0  =  q(qhq)  *  (49) 

By  applying  the  same  procedure  used  in  Section  A,  the  2-D 
DO  A  estimation  can  be  found  using  the  proposed 
orthogonal  projection  real  rank  revealing  LU  factorization. 
Note  that  the  difference  between  the  two  proposed  methods 
is  how  to  construct  the  null  spaces  to  estimate  the  azimuth 
and  elevation  angles  for  the  incident  sources 


III.  SIMULATION  RESULTS 

For  simulation,  the  spacing  between  the  two  adjacent 
elements  in  the  uniform  linear  array  was  set  to  a  half 
wavelength  of  the  incoming  signals.  Further,  one  single 
snapshot  per  trial  and  50  independent  trials  in  total  were 
tested.  The  performance  of  the  two  proposed  methods  (OP- 
RRRLU)  and  (OP-RRRQR)  were  tested  when  the  sources 
are  coherent. 

Figures  2  through  5  show  the  plots  of  2-D  DO  As  for 
the  two  proposed  methods  (OP-RRRQR)  and  (OP- 
RRRLU),  respectively.  We  considered  K=3  coherent 
sources,  SNR=70  dB,  and  eleven  elements  were  assumed 
in  each  antenna  array  (the  element  at  the  origin  is  common 
to  both  arrays).  The  elevation  and  azimuth  angles  of  the 
two  sources  were  set  to  (75°,  60°),  (90°,  80°),  and  (110°, 
95°)  for  sources  1,  source  2,  and  source  3  ,  respectively. 
Fifty  independent  trials  were  performed  for  each  figure. 
Figures  2  and  3  show  the  histogram  plot  for  the  elevation 
angle  and  the  pow  er  spectrum  plot  for  the  azimuth  angle, 
respectively,  for  the  proposed  (OP-RRRQR)  method.  It  is 
observed  clearly  in  Figures  2  and  3  that  the  proposed 
algorithm  gives  accurate  2-D  elevation  and  azimuth  DOA 
estimations  both  sources  and  an  exact  peak  occurred  at 
(75°,  60°),  (90°,  80°),  and  (1 10°,  95°)  in  almost  all  the  cases. 
Figures  4  and  5  show  the  histogram  for  the  elevation  angle 
and  the  power  spectrum  for  the  azimuth  angle  using  the 
proposed  (OP-RRRLU).  From  these  figures  we  observe 
that  (OP-RRRLU)  gives  an  accurate  estimation  for  the 
azimuth  and  elevation  angles  for  the  three  sources,  and  a 
clear  peak  occurs  at  the  exact  directions.  Note  that  the  two 


proposed  methods  can  estimate  the  2-D  DOA  estimation 
efficiently.  However,  the  complexity  of  (OP-RRRLLI)  is 
lower  than  that  of  (OP-RRRQR)  by  a  factor  of  tw  o. 


iy4^ONCLUSIONS 

In  this  paper,  we  propose  two  methods  for  estimating 
the  2-D  direction  of  arrival,  employing  the  Root  MUSIC 
and  MUSIC  algorithm  for  coherent  and  non-coherent 
sources.  The  received  data  is  arranged  into  a  Toeplitz 
matrix  which  enables  us  to  perform  the  estimation  using 
only  a  single  snapshot  and  detect  the  incident  sources 
whether  they  are  coherent  or  not  without  any  spatial 
smoothing.  In  addition,  the  unitary  transform  successfully 
employs  to  transfer  the  complex  data  of  RRQR  and  RRLU 
factorization  to  real  data.  Moreover,  it  does  not  require 
pair-matching  between  elevation  and  azimuth  angle 
estimation.  These  advantages  make  our  proposed  method 
suitable  for  real  time  implementation. 


Figure.  1.  Block  diagram  for  proposed  algorithm  for  joint  2-D  DOA 
elevation  and  azimuth  angle 


I  ^ 


RRRQR  Factorization  Method 

- 

with 

Unitary  Root  MUSIC 

- 

- 

o - i - , - > - i - - - J — , - , - 1 - 1 - i - 

0  15  30  45  60  75  90  105  120  135  150  165  180  195 

Elevation  (angle) 


Figure  6.  Histogram  of  elevation  estimations  in  the  2-D  DOA 
estimations  for  three  coherent  sources  at  [(75°,  60°),  (90°,  80°),  (110°, 
95°)]  using  the  (OPRRRQR)  Method 
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Figure  7.  Power  spectrum  of  azimuth  estimations  in  the  2-D  DOA 
estimations  for  three  coherent  sources  at  [(75°,  60°),  (90°,  80°),  (110°, 
95°)]  using  the  (OPRRRQR)  Method 


RRRLU  Factorization  Method 
with  Unitary  Root  MUSIC 


)  75  90  1  05  120  135 

Elevation  (angle) 


Figure  8.  Histogram  of  elevation  estimations  in  the  2-D  DOA 
estimations  for  three  coherent  sources  at  [(75°,  60°),  (90°,  80°),  (110°, 
95°)]  using  the  (OP_RRRLU)  Method 


Figure  9.  Power  spectrum  of  azimuth  estimations  in  the  2-D  DOA 
estimations  for  three  coherent  sources  at  [(75°,  60°),  (90°,  80°),  (110°, 
95°)]  using  the  (OP_RRRLU)  Method 
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